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Abstract 

We consider kink-antikink collisions in the TDSG system introduced by Speight 
and Ward in 1994. We find that the TDSG kink supports extra internal modes of 
vibration and this results in resonance effects of the kind seen for the continuum 4 
system. 
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1. Introduction 



Nonlinear systems which have soliton solutions have played an increasingly important role 
in physics. One such area in physics is condensed matter theory where the soliton has been 
used to model crystal discolations, charge density waves, magnetic domain walls, etc |J. 
In many cases it has been necessary to model interactions of the soliton with other objects 
such as impurities, defects, phonons, antisolitons, or other solitons. In (l+l)-dimensional 
field theories, these interactions have sometimes been quite spectacular, particularly if the 
1-dimensional soliton, or the kink, supports extra modes of vibration. Such a mode can 
be excited during an interaction and this leads to a "resonance effect" [[|. This effect has 
been seen in both continuous and discrete systems [1-8]. These resonance effects give rise 
to "windows" in the space of impact velocity, i.e. there are certain intervals of velocities in 
this space for which the kinks form a bound state (a bion) and others for which the kinks 
scatter off each other. 

This paper deals with kink-antikink collisions in the topological discrete sine-Gordon 
system (TDSG) introduced by Speight and Ward in 1994. This system is particularly 
interesting since it is a discrete system but with many of the features of the continuum sine- 
Gordon. In particular, it maintains the Bogomoln'yi bound. The resulting Bogomoln'yi 
equation is a first-order difference equation which can be easily solved. The solution is 
analogous to the well-known continuum sine-Gordon static kink solution. Moreover, as we 
show in section 4, the TDSG kink supports an internal shape mode. 

The paper is divided into 5 sections. Section 2 contains a brief summary of kink dy- 



namics in the TDSG system flIU |, [Tj| . In section 3 we consider the kink-antikink collisions 
in the TDSG. In section 4 we show that TDSG kink supports an internal mode of vibration 
and section 5 contains some concluding remarks. 



2. Kink dynamics in TDSG 

The TDSG system is defined by the Lagrangian density 
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where the kinetic energy is Ek = f X/Vn ano - the potential energy is £p = | XX-^n + ^n)- 
The quantities .D n and F n are given by 



2 1 1 

D n = - sin -(/0 n +i - V'n) and F n = sin -{ip n+ i + ip n ), (2) 



giving 



^ = 7 E w sin2 2^™ +1 ~ V"n) + sin2 + V>n). (3) 

■0 n = -^(x) at the n i/l lattice site and ip n +i = i[)(x + h), where h is the lattice spacing. The 
first term in equation (|3|) represents an attractive force analogous to the Hooke force. The 
second term is the substrate potential, found by taking an average of the potential of the 
two nearest neighbours. 

In the continuum limit, as h — > 0, D n —> ip x , and F n — > sin^, and one recovers the 
standard expression of the continuum potential 



1 r°° 

Ep*^) = -7l (^ + sin 2 ^)dx. (4) 

^ J — oo 

The key feature of this model is the choice of the derivative on the lattice. This choice is 
motivated by the Bogomol'nyi argument of the continuum system. Rather than choose the 
standard forward difference, Speight and Ward choose, 



h 

with the kink boundary conditions 



COS^ n+ i - COSljjn 

= D n F n , (5) 



lim ip n = and lim ip n = it. 



n— >— oo 
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This leads to the following discrete Bogomol'nyi argument 



= E P + ^ J^(cos ip n+1 - cos ip n ) (7) 

= E P - 1, (8) 
so E P > 1, with equality if D n = F n , i.e. if 
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-sin-(/0n+l - Ipn) = Sin-^n+l + V'n) (9) 



or, 



■0 n = 2tan 1 exp a(n/i — 6). (10) 



The constant a = /i" 1 ln[|^|] is the kink's slope and b G R is the position. Equation ( |10|) is 
the static kink solution. So the moduli space of static solutions for the TDSG is isomorphic 
to R rather than Z, which is what one usually expects for a discrete system. Moreover, Ep 
is independent of b suggesting that the model has no Peierls-Nabarro barrier. This means 
that the kink can move along the lattice arbitrarily slowly without getting pinned. 
The Euler-Lagrange equation gives the following equation of motion, 



4 - h 2 4 + h 2 
= Ah 2 cos ^'"( sin ^"+ 1 + sin^ n _i) ^-sin^ n (cos?/' n+ i + cos^ n _i). (11) 

Equation (|TTD can then be used to study the behaviour of a moving kink. Speight and 
Ward have solved fliTD using a 4th-order Runge-Kutta scheme. Their initial condition was 
a static kink Galilean boosted to a velocity v. Their results show that for small initial 
velocities the kink wobbles with a period h/v as it moves along the lattice. This is a purely 
dynamical effect since there is no Peierls-Nabarro barrier. 

For fast moving kinks it is found that the kink radiates as it moves along the lattice. 
This is observed through a gradual decrease in the kink velocity. But it is found that 
there is certain velocity threshold below which the radiation from the kink is significantly 
reduced. 
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The existence of such a velocity threshold can be understood by considering the lin- 
earised equation of motion |L0| . 



4 - h 2 4 + h 2 

Equation (|l"2D can be used to derive a dispersion relation for small amplitude travelling 
waves, given by 



4 + h 2 4 - h 2 

2^ C ° S ^' (13) 

where k is the wavenumber. Since < h < 2, we have 1 < u < 2/h. The case uo = 1 
corresponds to the threshold velocity. The frequency at which the kink hits the lattice sites 
is v/h per unit time. So provided v > h/2-K, the kink will radiate (at uo = 2irv/h). But if 
v < h/2n then u < 1 and the kink does not radiate. Hence there is a "preferred" velocity 
for the kink. Since this velocity depends only on h, it can be regarded as a feature of the 
lattice. 

The choice of the functions D n and F n is not unique. For instance, D n could be 
multiplied by a function f(x, h) and F n by f(x, h)^ 1 and the product is independent of / 
so one would still have a topologically stable kink solution. The only condition on f(x, h) 
is that lim/j^o / — 1- Similarly there is freedom in the choice of ip. Zakrzewski has 
investigated the effects on the kink motion when this choice is exploited. 

The ip in (|ll|) is replaced by 



i) -> ^ (14) 

1 + cr?/> 2 



where a = a(h) and a(0) = 1. The equation of motion QlTD is now changed to 



■ 2 

(1 + aijj n ) 



P 
h 2 



(sin^n+i - ^n) ~ sin(^n - V'n-l)) 



4/ 2 



(sin(^n+i + ipn) + sm(ip n + ip n -i))- 

(15) 
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The choices of / and a are motivated by the fact that equation QT3 ) still needs to admit 
the Galilean boosted sine-Gordon kink solution. One choice is the following 



1-1)2 



,2 




) 



2 



+ 



4sinh 2 

— W 



M. 



(16) 



8 sinh 2 



where /? = /i/ 2 \A 



— vf, and 



a(h) 




sinh 2 (2/3). 



(17) 



Numerical simulations of (|T5|) using a 4th-order Runge-Kutta method has shown that 
the kink propagates without emitting any radiation if the initial velocity of the kink is 
chosen to be the velocity v\ used in the expression of f 2 . For initial velocities v < t>i, there 
is virtually no emission of radiation so the kink velocity hardly decreases. For v > V\, the 
kink initially propagates at velocity greater than v± but then slows down to a velocity just 
less than V\ after which it's virtually constant with little emission of radiation. The initial 
increase in velocity is due to the readjustment of the kink configuration to match the lattice 
field configuration. 

3. Kink-antikink Collisions in the TDSG System 
3.1 Preliminaries 

In this section we present results of kink-antikink interactions in the Speight-Ward model. 
The results are obtained by solving (|TTD with the initial condition corresponding to Galilean 
boosted static kink and antikink solutions at each end of the lattice grid, i.e. 



ip(x, 0) = 2 tan 1 exp a(x + b) — 2 tan 1 exp a{x — b) 
ip(x, 0) = — av(sech a(x + b) + sech a(x — &)). 



(18) 
(19) 
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The boundary condition is taken to be ip(x m i n ) = ip(x max ) = 0, i.e. the ends of the kink 
and anti-kink are held fixed at all times. 

The kink and antikink move towards each other with relative velocity 2v. Equation (|TT|) 
is solved using a 4th-order Runge-Kutta algorithm. The programme is run with various 
lattice spacings h. The time step is chosen to be 0.05. This conserves the energy to within 
0.05% in a simulation which runs for 1000 units. We use the same time step for all values 



Our choice of the boundary condition means that any radiation emitted by the kink 
will reflect off the edges and eventually interact with the kinks. This is easily avoided by 
making the grid sufficiently large. 

As mentioned already, the amount of radiation emitted by the kink depends on the 
kink velocity and also on the lattice spacing. For large velocities and coarser lattices the 
amount of radiation is significantly greater. In these cases we have damped the first few 
lattice sites at each end of the grid. That is, the value of ip at these sites is decreased by a 
constant amount (10%) for the entire duration of the simulation. The total energy of the 
system is of course no longer conserved. The energy loss depends on the initial velocity of 
the kink and on the lattice spacing, but it is usually within 10% of the initial energy, for a 
simulation running for 1000 units. 

3.2 Simulation Results 

The simulations showed that the outcome of an interaction depends on the impact velocity 
of the kinks. The velocity of the kink (or the antikink) is defined to be the rate of change 
of X(t), where X(t) is the "average" position of the kink (or the antikink). X(t) for the 
kink is given by 



of h. 



o 




X(t) 



i=-N 



(20) 



o 




i=-N 
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Here ej is the total energy of the field and the radiation at the i lattice site, where 
Xi = x m i n + ih. The quantity Y1°-n e « * s the total energy of the field and the radiation 
for the single kink. We found that for cases when the phonon radiation from the kinks 
is not significant, X(t) makes a good approximation to the kink position. We define the 
velocity of the kink to be X(t), which is computed using the forward difference, X(t) = 
±(X(t + At)-X(t)). 

Since the kink wobbles as it moves along the lattice, X(t) is an oscillating function. The 
period of oscillation is the wobble period, h/v, where v is the velocity used in the initial 
condition ([19j). This wobble can be understood by considering the geodesic motion on the 
kink moduli space. This moduli space has a periodic metric [|I(J . 

Figure 1 shows various plots of kink velocity against time for kinks with different initial 
positions. The kinks are on a lattice with h = 1.4 and the initial velocity of the kink is 
v = 0.197. In all cases the average velocity of the kink is significantly decreased. This is 
due to phonon radiation from the kinks. The plots also show that the kink velocity depends 
on the initial position of the kink, modulo nh, where n € Z. 

The outcome of an interaction has sensitive dependence on the velocity of the kink. In 
particular, there are three different outcomes: 

(a) If the kink impact velocity is above a certain critical velocity, v c , then the kinks 
simply pass through each other, i.e. there is a smallest number v c such that the outcome 
of the interaction for all v > v c is always a passing through behaviour. The value of v c 
depends on the lattice spacing and also on b. Figure 2 shows the dependence of v c against 
h. The dotted curve represents the quadratic 0.073/i 2 + 0.036/i. This curve is an empirical 
fit to the data points and is not derived from the theory. 

The dependence of v c on the position of the kinks is due to dynamical dressing and the 
kink wobble. For small initial velocities (figure 1), the kink velocity oscillates with a period 
{h/v). The amplitude of oscillation is of the order of 10~ 2 . For large velocities, the phonon 
radiation make the oscillations erratic and also significantly decrease the overall velocity of 
the kink. The initial velocity also changes due to the readjustment of the kink field to suit 
the lattice distribution (dynamical dressing). The initial propagating velocity of the kink 
is not always v, the velocity used in expression (0). In some cases it is larger than v and 
in others it is less than v. This change depends on the magnitude of v and on the lattice 
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Figure 1: A plot of the kink velocity against time. The three curves represent velocities of 
kinks with an initial velocity of 0.197 but starting at different points on the lattice. The 
solid curve represents a kink with b = 30, the dashed curve a kink with b = 25 and the 
dotted curve is for b = 15. 



spacing. The effect can be seen in figure 1. 

For initial kink velocities which give impact velocities below v c , we found that the 
outcome is either (b) a bounce or (c) a trapped breather-like state (a bion). The bounce is 
essentially a kink-antikink scattering, where the kinks after the first collision pass through 
each other but fail to escape to infinity. Instead, they travel a small distance and then turn 
around for another collision. After the second collision the kinks either escape to infinity 
or turn around again for more collisions. For a 2-bounce event (meaning the kinks pass 
each other twice) the ultimate outcome of this scattering would be a reflection, and for a 
3-bounce a transmission. 

The velocities for the bounces occur in small intervals. For velocities outside of these 
intervals the outcome is always a bound state. This is similar to the results for the </> 4 model 
f|, except in this case the outcome does not depend fractally on the impact velocity. We 
have only seen 2, 3 and 4-bounce events. But the velocities at which the bounces occur are 
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Figure 2: A plot of critical velocity against lattice spacing. The crosses represent data 
points found numerically. The curve represents a theoretical fit of v c = 0.073/i 2 + 0.036/i. 



not easy to find and the intervals are small (10~ 7 for a 4-bounce), so it may be that higher 
bounce events occur but we have not seen them. 

Figure 3 shows an example of the 2-bounce event. In this figure we have plotted the field 
ip(0,t) against t. ip(0,t) is the field at the centre of mass of the kink-antikink system.The 
kinks are on a lattice of unit spacing with an initial separation of 10 units. v c for this 
combination is found to be 0.100. The kinks bounce for all initial velocities in the range 
0.0969 < v < 0.0972 (2-bounce window). 

Figures 4 and 5 show examples of three and four bounce events. In figure 4 the kinks 
transmit whereas in 5 they reflect off each other. As can be seen from the figure for the 
2-bounce case, the kink and antikink collide, pass through each other and travel a small 
distance before stopping, turning around and travelling in the other direction. The kinks 
then collide again, pass through each other and this time escape to infinity. Figure 6 is a 
graph of ^(0, t) against t representing the oscillations at the centre of mass of a bion. This 
is essentially an n-bounce event for n very large. The erratic oscillations suggest that the 
system might be chaotic. The period of oscillations become constant after about 10 cycles 
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Figure 3: A figure representing the 2-bounce event. The kinks collide and essentially reflect 
off each other. 
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Figure 4: A 3-bounce event where the scattering of the kinks results in a transmission. 
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Figure 5: A plot showing a 4-bounce event. 
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Figure 6: A plot of ip(0,t) against time, representing the oscillations at the centre of mass 
of the bion. 
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and the bion then oscillates for a long time (i.e for hundreds of cycles). We have calculated 
the maximum Lyapunov exponent of the time series -0(0, t), using the algorithm of Wolf et 



al |L5|. Its value is approximately 0.18, suggesting that the oscillations are mildly chaotic. 
In general, a positive maximum Lyapunov exponent is taken as a formal definition of chaos. 

The results of the simulation can be understood using the resonance energy exchange 
mechanism. This requires the TDSG kink to support shape modes. We show next that 
the TDSG kink does support such a mode. The frequency of this mode is calculated in 
two ways, by using a collective coordinate approximation and also numerically, by Fourier 
analysing the deformations of a perturbed static kink. Note that we could also have used 



a matrix-based approach similar to the one used by Speight in ||16|| , where it was used to 
calculate the Casimir energy of the kink. 

4. Kink internal shape mode 

4.1 Collective coordinate approximation 

We can obtain a collective coordinate approximation of the TDSG system by treating the 
scale factor a in the kink as a dynamical variable. This reduces the number of degrees of 
freedom from infinity to one. The Lagrangian of the reduced system is given by 

C = f( ai b)d 2 -g(a,b), (21) 
where f(a, b) and g(a, b) are functions given by 

/(a, 6) = ^J2(x-b) 2 sedi 2 a(x-b), (22) 

X 

g(a,b) = i tanh(^a/i) + ^coth(^aZi). (23) 
The corresponding equation of motion is 

2/fi + fa 2 -g' = 0. (24) 
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where the prime denotes differentiation with respect to a. g(a, b) has a minimum (stable 
equilibrium) at a = s := |tanh _1 |. Small amplitude oscillations a(t) = s + e(t) satisfy 
the linearised equation, 

2f(s,b)e = g"(s,b)e. 

Hence the waves oscillate with frequency 

= J_ / g"(s) 
V 2n\j2f(s,by 

g"(s) = (1 — \h 2 ) 2 , so the frequency of the shape mode depends on the lattice spacing. For 
h = 1, the shape mode frequency is 0.152. Table 1 shows the dependence of frequency on 
the lattice spacing h. 



h 


frequency 


0.1 


0.175 


0.2 


0.175 


0.5 


0.170 


0.8 


0.161 


0.9 


0.157 


1.0 


0.152 


1.2 


0.143 



As a check on the accuracy of the collective coordinate approach, we have used it to 
calculate the frequency of the shape mode of the 4 kink. The value found is ^/(1.55), 
which is fairly close to the actual value a/(1.5). 

4.2 Numerical method 

Numerically, the shape mode frequency is obtained by Fourier analysing the deformations 
of the perturbed static kink. We set the initial condition for equation ([□]) to be 



(25) 
(26) 



if>(x,0) =tp K {x,0) + e 
tP(x,0) = 0, 
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(27) 



where iI>k(x) is the static kink solution and e is the pertubation. The value of e is chosen 
to be -0.03. The field ip(x,t) is sampled at various values of x. From this we obtain the 
variations in the field, 5ip(x,t) = ip(x,t) — ip(x,0), for each of the sampled values of x. 
5ip(x,t) were then Fourier analysed using the MATLAB FFT routine. 

In figure 7 we have plotted the power spectrum for the data sampled at x = for a kink 
on an h = 1 lattice. The power spectrum of the data for other values of x give identical 
power spectrums. 
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Figure 7: Power spectrum of the fluctuations around a perturbed static kink 



Figure 7 was constructed using N = 2 14 = 16384 data points. Aliasing of the power 
does not occur since the amplitude at the Nyquist frequency is essentially zero. There is a 
peak at v — corresponding to the Goldstone mode, and a peak at v — 0.161 corresponding 
to the shape mode. 

This procedure is repeated for other values of h. In contrast to the results suggested by 
the collective coordinate approach, we find that the frequency of the shape mode is more 
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or less independent of h. For 0.5 < h < 1.2, we found 0.158 < f < 0.161. The value of the 
frequency (for h — 1) however is within the range of frequencies found using the collective 
coordinate method. 

In the continuum limit the TDSG shape mode frequency is precisely the quasimode 



frequency of the continuum sG system found by Rice, ujr = 0.175 [12|, although it is 
unlikely that such a quasimode exists [ITT 



It is interesting to compare the situation with other discrete sine-Gordon systems. The 
well-known Frenkel-Kontorrova model does not support a kink with an internal shape mode 
for all values of the discreteness parameters. The sine-Lattice model (s-L) however does 
support one, as long as the system is only weakly discrete |7[] . The shape mode for the s-L 
system lies just above the lower phonon band suggesting that it is a genuine kink shape 
mode rather than just a phonon resonance. Moreover, in the continuum limit this shape 
mode converges to the sine-Gordon quasimode frequency. The existence of a shape mode 
in a discrete sine-Gordon system depends on the anharmonicity of the potential term. 



5. Conclusion 

We have found that the kink-antikink interactions in the TDSG system exhibit resonance 
phenomena similar to those in the continuum 4 system. These are found to be due to the 
excitation of an intermal mode of the TDSG kink. A collective coordinate analysis shows 
that the frequency of this mode depends on the lattice spacing, and in the continuum limit 
{h — > 0) the frequency seems to approach the quasimode frequency predicted by Rice. 
The recent analysis done in ]17| however says that there cannot be a quasimode for the 



continuum sine-Gordon kink, but they do predict the existence of a quasimode for a lattice 
sine-Gordon kink. This is not the quasimode we have found, so this is something genuinely 
new. Further numerical simulations are required to determine what is happening in the 
limit h — > 0. 



Also, as we have mentioned already, there is a modified TDSG system |]I4[ which admits 
an exact travelling-kink solution (with a fixed velocity). It would be interesting to perform 
kink-antikink interactions for this system. Since the expressions we have used to perform 
the interactions are only an approximation to the equation of motion, it would be interesting 
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to see the difference. 

We have also done kink-antikink collision simulations for the TD</> 4 system ||18|| . The 
results are found to be analoguos to the conventional discrete 4 system [|J, though in the 
"windows" structure in the space of impact velocity is not a fractal. 
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